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rS^ • Abstract 

We have measured the inclusive D r± production cross section in a two-photon collision at the 
TRISTAN e + e~ collider. The mean \fs of the collider was 57.16 GeV and the integrated luminosity 
was 150 pb' 1 . The differential cross section (da(D* )/dPr) was obtained in the Pt range between 1.6 
and 6.6 GeV and compared with theoretical predictions, such as those involving direct and resolved 
photon processes. 

Hadron production in two-photon collisions is described by the vector-meson dominance model, the 
quark-parton model (direct process) Rb and the hard scattering of the hadronic constituents of almost-real 
photons (resolved photon process) j^Tg, [§ ||, which has been observed by the previous experiments 0. 
However, more detailed studies are necessary in order to understand these processes quantitatively. Heavy 
quark pair production processes are good probes, since the theoretical calculations are less ambiguous 
than for light quarks 0. 

The previous measurements had been carried out at around yfs ~ 30 GeV || [|, [l(|, and are consistent 
with a recent theoretical prediction |7J]. At the TRISTAN energy (^/s ~ 60 GeV), the cc production cross 
section becomes sizable; we have obtained the largest statistics for this type of process. We carried out a 
measurement of the D*^ production cross section at a Pt greater than 1.6 GeV in two-photon collision 
events. 

The detail concerning the TOPAZ detector can be found elsewhere |ll[. The integrated luminosity 
of the event sample used in the analysis was 150 pb~ l . The mean < y/s > of the collider was 57.16 GeV. 
The trigger conditions are as follows: more than two tracks with Pt > 0.3 ~ 0.7 GeV and opening angle 
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> 45~70 degrees (depending on the beam condition); a neutral energy deposit in the barrel calorimeter 
be greater than 4 GeV; or that in the endcap calorimeters be greater than 10 GeV. 

The selection criteria for two-photon events are as follows: the number of charged tracks be > 4; the 
total visible energy in the central part of the detector be between 4 and 25 GeV; the vector sum of the 
transverse momenta of the particles with respect to the beam axis (|EPt|) be less than 7.5 GeV; the 
sum of the charges be < 3; the event vertex be consistent with the beam crossing point; and no large 
energy clusters (E > 0.25Ei, eam ) in the barrel calorimeter. In addition, we divided each event into two 
jets with respect to the plane perpendicular to the thrust axis, and the cosine of the angle between the 
two jets was required to be greater than -0.9. These restrictions were made in order to reduce beam-gas, 
and single-photon-exchange events. A total of 10788 events was selected. 

The charged track selections were as follows: the closest approach to the event vertex be consistent 
within the measurement error; the number of degrees of freedom (DOF) in the track-fitting be > 3; and 
P T be > 0.15 GeV. 

The 7 selections were: the cluster be detected by a barrel-type lead-glass calorimeter; the energy be 

> 0.2 GeV; and the cluster position be separated from any charged-track extrapolations. In addition, the 
e + e~ pairs which were consistent with 7 conversion at the inner vessel of the Time Projection Chamber 
(TPC) were reconstructed (lC-fit) and used as 7 candidates. 

In order to reconstruct charm quarks, we used the decay modes D*^ — > n ± D°(D ), followed by 
(D°(D°) — » K^ir X(X)). From now on, any mention of decays includes the charge conjugation modes, 
for simplicity. For D° decays, the decay modes D° — ► K~n + , K~ir + p°, K~n + ir° (K*~n + , K~p°, 
K*°ir°), and K~ir + ir°ir° were reconstructed by using kinematical constraint fits (1-3C). The cuts on x%t 
were required to be greater than the 5% confidence level (CL). These decay modes were selected because 
they have relatively high acceptances, considering the branching fractions and detector acceptances. For 
more than two-body decay modes, we applied a dE/dx cut in selecting K~; for a two-body decay, 
a sufficient S/N was obtained without this cut. The cuts on the vector meson masses were carried 
out according to the detector resolution and their intrinsic decay widths. In the case of vector-plus- 
pseudo-scalar decay, we applied cut on the angle of the decay product of the vector meson (8vp) in its 
center-of-mass frame with respect to the vector meson line of flight, i.e., |cos#yp| > 0.5. 

The D* s were reconstructed with those D° candidates mentioned above while combining 7r + s (soft- 
pions from hereafter) with momenta less than 0.65 GeV. The energy fraction, z = E(D* + )/Eb ea ,m, was 
required to be between 0.1 and 0.25. We then calculated the mass differences, i.e., dM — M(ir + D°) — 
M(D°). The dM distribution is plotted in Fig. 1(a). 

There were multi-D°-candidates for one soft-pion which gave similar dMs. These occurred when one 
of the lowest momentum particles of the D° decay was misidentified|13]. In these cases, each entry gave 
similar dM value. The differences of dM values for pairs sharing the same soft-pions in the same events 
were plotted in Fig. 0. The peak around zero caused an overestimate of the statistical significance of 
the signal. In order to cure them, we carried out a weighting method, i.e., when there was more than 
one D°-candidate for a given soft-pion, each entry in the dM histogram was weighted by the reciprocal 
of the number of candidates. By this procedure, we could avoid any overestimation of the statistical 
significance of the peak entry. This was checked using a Monte-Carlo simulation. The resulting dM 
distribution is shown in Fig. Q-(b). The mean visible energy at the rest frame (WVIS) for the events 
containing the D*^ candidates were 5.3 GeV. In order to check the peak around the D*^ region, we 
carried out the wrong-sign combination such as D°tt~ and etc. The results are plotted in Fig. [j]-(c), (d), 
and (e). There were no peak structures at all. The excess below the D* ± peak was explained by the 
energy loss at support structures of the inner field cage and the central membrane of the TPC, which were 
distributed inhomogeneously, whereas the correction for them were made only in average. Figs |3| show the 
mass-differences of two cases: (a) when the soft-pion passed near these supports, and (b) when it passed 
far away from them. This happened because the soft-pions in the two photon events had extremely low 
momenta. We counted the entries of the higher dM peak as the D yield in the experiment, and in 
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Figure 1: Mass-difference (dM = M(£)°7r ± ) - M(D )) distributions: (a) dM distribution without the 
weighting; (b) after the weighting. The curve is obtained from a best-fit function described in the text; (c) 
dM distribution for wrong-sign soft-pion; (d) that for wrong-sign kaon; (c) that for wrong-sign soft-pion 
and wrong-sign kaon; and (f) background estimation from e + e~ — > 7 — > qq processes. 
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Figure 2: Differences of dM values for pairs sharing the same soft-pions in the same events. 



the Monte-Carlo we corrected the peak entries about -15% considering the amount of the materials and 
their solid angle from the interaction point. The error of this correction was considered to be 5% which 
was estimated from the e + e~ conversion pair yield in the two photon sample. Concerning the higher dM 
fluctuation, we have no explanations other than statistical fluctuations. The fitting function was the sum 
of a Gaussian and the following function: 

a(dM - M(tt ± )) 6 (1 - dM/c) n , 

where a, b, and c are free parameters and n is an average multiplicity of the event sample. The obtained 
D* ± yield was 20.0 ± 5.0, where the x 2 of the fit was 31 with DOF=35. The peak position and its width 
obtained by this fitting were 147.4 ± 0.2, and 0.8 ± 0.2 MeV, where the detector simulation predicted 
145.4 and 1.6 MeV, respectively. We considered that the shift of the peak position was caused by the 
inhomogeneous material distribution described above, whereas the energy loss correction was carried 
out assuming uniform material. The shift of the mass-difference peak quantitatively agreed with the 
expectation. The resolution difference was due to the overestimation of the material in front of the TPC 
in the simulation in the most probable energy-loss case. 85% of the soft-pions were considered to pass 
away from the support structures of the TPC. We tried x 2 and likelihood fitting and also tried polynomial 
background functions. The differences in the total D*^ yields were within 10%. Thus we concluded that 
the systematic error of the fitting procedure was 10%. The yields for the four decay modes described 
before were 9.5 ± 3.3, 4.8 ± 2.3, 4.8 ± 2.2, and 0.9 ± 0.9, respectively. 

The background from e + e~ — > 7 — > 99(7) distributed smoothly in this dM plot (height was about 1 
count a bin as shown in Fig. Q-(f)). No peak structures were obtained by the Monte-Carlo simulation. 
We also checked the event vertex distribution in order to determine the contamination of the beam-gas 
background. There was a clear peak at the interaction point with no tails. 

In order to compare the experimental data with the theoretical prediction, we chose da(D* )/dPr, 
instead of the total cross section, as had been reported |g, g, [t0|, since the accepted D*^ events were 
limited to the high Pt region due to the detector acceptance. We were sensitive to those D* ± s with 
transverse momenta (Pt) between 1.6 and 6.6 GeV. The lower limit was due to the detector acceptance 
and the higher due to statistics. An acceptance correction was carried out by using a lowest-order (Born 
approximation) direct-process Monte-Carlo simulation, in which an equivalent photon approximation was 
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Figure 3: Mass-differences of two cases: (a) when the soft-pion passed near the TPC support structures, 
and (b) when it passed far away from them. 
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Figure 4: da(e" 



e D X)/dPr versus Pt- The data points with error bars are the experimental 



data and the curves were obtained by Monte-Carlo simulations while assuming various processes. The 
hatched area is a prediction based on the direct process, the solid curve is based on a combination of 
direct and resolved photon process by the GRV parametrization, the dotted one is based on the LAC1 
parameterization. The vertical scales for the resolved photon process were normalized by the relative 
acceptance to the direct process, because of the acceptance difference described in the text. 



used. The Pt(7?* ± ) smearing by the Pt of the initial photons was checked by a full calculation of the 
e + e~ — > e + e~77 — > e + e~cc process, and was concluded to be small. The hadronization was taken care 
of by LUND6.3 Monte-Carlo H M- We did not use a parton shower option in the LUND Monte-Carlo, 
since the choice of Q 2 of the event was not well defined. The fragmentation parameters in the calculation 
have been tuned to reproduce the general properties of single-photon exchange processes. The systematic 
ambiguity for 73* ± inclusive spectrum was considered to be at most 10 %. Because the fragmentation of 
the D was measured well at various y/s and the symmetric LUND fragmentation function reproduces 
it very well. The systematic uncertainty of the cross section is mostly due to the fitting procedure, 
fragmentation scheme (described so far), the branching ratios of 73° |1J] and D* + |1(|, the ambiguities 
in the detector acceptance and the calibrations. It was estimated to be 15%. The resulting differential 
cross sections are 57 ± 24, 9.5 ± 5.1, 6.3 ± 3.2 and 2.0 ± 1.8 pb/GeV (systematic errors included) for Pt 
regions 1.6-2.6 ,2.6-3.6, 3.6-4.6, and 4.6-6.6 GeV, respectively, and are plotted in Fig. [|, together with 
the theoretical predictions. The raw entries at each Pt binnings were 7.2±2.9, 6.4±3.3, 4.0±2.0, and 
2.4±2.1, where the errors are statistical error only, respectively. 

Theoretical predictions for both direct and resolved photon processes were compared with this mea- 
surement as follows. In order to estimate a(D* ± ) from the theoretical calculation of a cS , we assumed 
P C ^ D ,± x B(D*+ -> n+D ) = 0.185 ± 0.024 at ^5 = 90 GeV in Ref. @. This value is consistent 
with that obtained at y 7 ! = 10 GeV Q. The events were generated by the lowest-order calculations. 
The relative acceptance for the process of the resolved photon to that of the direct process was obtained 
to be 0.76 ± 0.25 due to the event selection, where the error was especially due to the gluon-density 
function dependence. Although, this ambiguity is small compared with the total contribution including 
the direct process. Since we did not measure the energy flow in the low-angle region, spectator jets were 
not detected. This lower acceptance was especially due to the lower visible energy cut of 4 GeV. The 
vertical scale for the cross sections of the resolved photon process shown in Fig. was normalized by this 



factor. 

Corrections of order a s were carried out according to the procedure in Ref. JlOJ ] . Since string frag- 
mentation includes a parton-shower-like effect, the next-to-leading-order correction in the Pt spectrum is 
doubly counted. We used instead the Pj^-independent correction to the direct process. The cross section 
of the direct process is increased by a factor of 1.31 uniformly and that of the resolved photon process 
is corrected by a Pt dependent function, due to the presence of the process 7(7 — ► ccq (a part of this is 
absorbed in the gluon density function in the resolved photon) . The Pt dependent factors were obtained 
by the following way: At first we derived the Pt dependent ratios between the higher and the lowest 
order calculations for the direct and resolved photon processes; we then calculated the factors between 
those of the direct and resolved photon process; and those factors were normalized to fit to the total 
cross section of the higher order calculations for the resolved photon process. The resulted Pt dependent 
correction were written as 

0.50P T + 0.54, 

where P T is a Pt of charm quark. The charm quark mass (m c ) and the renormalization scale (/x) 
were assumed to be m c =1.6 GeV and \i = \/2m c , respectively. The curves in Fig. || represent these 
predictions. The parametrization dependence for the resolved photon process appears in the lowest Pt 
binning. Although our data favor GRV|| or LAClQ, none of them explain the high Pt excess (+1.5ct 
at P T > 2.6 GeV). 

There are ambiguities in the above-mentioned theoretical predictions, which depend on m c and fi. We 
changed m c from 1.3 to 1.8 GeV and fi from m c /2 to 2m c . The lower m c and fi give a higher cross section 
which has a Pt dependence. The lowest Pt region may be explained by such ambiguities as mentioned 
above, and the threshold enhancement of cc production. We thus concentrate on the higher Pt regions, 
i.e., Pt >2.6 GeV. The case for m c — 1.3 GeV and /i = m c /2 gives the highest one, i.e., +0.54pb (+6%) 
higher than the nominal value. The experiment gives <r(2.6 < Pt < 6.6 GeV) — 19.8 ± 7.0pb, where the 
GRV parametrization predicts 9.0 ± 0.5 pb. Then the difference is 10.8 ± 7.5pb. 

In addition, if gluons inside the photon have Pts of 0.4 GeV in average, the PtS of ccs are shifted. 
This may increase the resolved photon cross section by about a factor of two, because the spectrum of 
this process is proportional to an exponential function. Therefore, the 2 ~ 3pb increase may be expected 
by this (no increase in the direct process). Then our results become consistent with the expectation 
within 1 a level. 

Such event shapes as missing Pt and thrust distributions were checked. These shapes are consis- 
tent with the prediction of the lowest-order direct and resolved photon processes, except for the overall 
normalization factor. 

In a part of the data set (integrated luminosity of 90 pb^ 1 ), there were forward calorimeters (FCL: 
made of BGO) which covered the polar angle region between 3.2 and 13.6 degrees |2(J. We analyzed 
tagged events by the FCL with the same analysis. We observed the total D*^ yield to be 3.5 ± 3.6, where 
the lowest-order direct process predicts 1.0 ± 0.4 events. We therefore need more statistics in order to 
observe tags from electrons and positrons as well as spectator jets. 

In conclusion, we measured the inclusive D production cross section in two-photon collision events 
at the TRISTAN e + e~ collider. The mean y/s was 57.16 GeV and the integrated luminosity was 150 
pb^ 1 . We observed 20 ± 5 D*^ in our data sample. The production cross section in the high Pt region 
was <7d«±(2.6 < Pt < 6.6 GeV) = 19.8±7.0pb. Comparisons with the theoretical prediction of the direct 
process and the resolved photon process were carried out. 
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operation of TRISTAN. We also thank all engineers and technicians at KEK as well as members of the 
collaborating institutions: H. Inoue, N. Kimura, K. Shiino, M. Tanaka, K. Tsukada, N. Ujiie, and H. 
Yamaoka. 
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